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ABSTRACT 



Dusty hyper luminous galaxies in the early universe provide unique environ- 
ments for studying the role of massive stars in the formation and destruction of 
dust. At redshifts above ~ 6, when the universe was less than ~ 1 Gyr old, dust 
could have only condensed in the explosive ejecta of Type II supernovae (SNe), 
since most of the progenitors of the AGB stars, the major alternative source of 
interstellar dust, did not have time to evolve off the main sequence. In this paper 
we present analytical models for the evolution of the gas, dust, and metals in 
high redshift galaxies, with a special application to SDSS Jl 148+5251 (hereafter 
J1148+5251), a hyperluminous quasar at z = 6.4. Ignoring accretion of inter- 
stellar matter onto grains, the main free model parameters are the dust yield in 
SNe, and the grain destruction efficiency by supernova remnants. We find that 
an average supernova must condense at least 1 Mq of dust to account for the 
observed dust mass in J1148+5251. Theoretically, this large yield can only be 
attainable if stars are formed with a top heavy initial mass function. Observa- 
tionally, it is in excess of the largest dust yield of < 0.02 Mq found thus far in 
the ejecta of any SN. If future observations find this to be a typical supernova 
dust yield, then additional processes, such as accretion onto preexisting grains, 
or condensation around the AGN will need to be invoked to account for the large 
amount of dust in this and similar objects. The galaxy's star formation history 
is still uncertain, and current observations of the gas, metal, and dust contents 
of J1148+5251 can be reproduced by either an intensive and short burst of star 
formation {ip > 10^ Mq yr~^) with a duration of < 10® yr, or a much lower 
star formation rate {ip ~ 100 Mq yr~^) occurring over the lifetime of the galaxy. 
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Analysis of the spectral energy distribution of J1148+5251 suggests an AGN lu- 
minosity of ~ 7 X 10^^ Lq, requiring the early formation of a supermassive black 
hole of ~ 10^ Mq. 

Subject headings: galaxies: formation, evolution, high-redshift, starburst, AGN 
- infrared: galaxies, general - 

ISM: interstellar dust - individual (SDSS J114816.64+525150.3) 



INTRODUCTION 



Quasars and galaxies d etected at redshifts z > 6 and observed at f ar infrared (IR) and 



submillimeter wav elengths (iHughes et al.l Il998l : iBertoldi et al.l l2003al : iRobson et al.l |2004J : 
Beelen et al.ll2006l ) exhibit luminosities in excess of ~ 10^^ Lq, and inferred d ust masses and 
star formation rates in excess of ~ 10® Mq and of ~ 10^ Mq yr~^, respectively (IBertoldi et al. 



2003al ). For comparison, the Milky Way contains only about 5 x 10'' Mq of dust, for an as- 
sumed gas mass of 5 x 10^ Mq and a dust-to-gas mass ratio of ~ 0.01. Some of these high-^; 
galaxies are younger than ~ 1 Gyr, strongly suggesting that the dust must have condensed 
in core collapse supernovae which inject their nucleosynthetic products and newly-condensed 
dust into the interstellar medium (ISM) relatively promptly after their formation. In con- 
trast, AGB stars, the other main sources of interstellar dust, inject their dust products 
into the ISM only after their progenitor stars have evol ved off the main sequence, result- 
ing in a > 500 Myr delay time since their format ion (IDwekl Il998f). This delayed injec- 
tion o f AGB-condensed dust has been suggested by iDwek et al.l (120051 ) and iGalliano et al. 



(120071 ) as a possible explanation for the observed trend of the IR emission from polycyclic 



aromatic hydrocarb on (PAH) molecules with galaxies' metallicity (lEngelbracht et al.l 12005 
Madden et aUbood ). 



Supernovae (SNe) play a dual role in the evolution of interstellar dust. On the one 
hand, they are potentially the most important source of interstellar dust. If all the refrac- 
tory elements precipit ate with a 100% efficiency from the gas, a typic al 25 Mq SN can form 
about 1 Mq of dust (IWoosley fc Weaverl Il995l : iNomoto et al.l 120061 ). On the other hand, 
their expanding blast waves sweep up interstellar dust grains, destroying thena by thermal 
sputtering and evaporative grain-grain collisions (I Jones et al.lll996l : |Jonesll2004l ). The pres- 
ence of the large quantities of interstellar dust in these youthful galaxies therefore provides 
a unique opportunity for studying the dual role of SNe in the formation and destruction of 
dust. 



The importance of SN condensates as important contributors to the abundance of inter- 
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stellar dust was previously recognized by many authors ( 


Dwek & Scalolll980 


Eales &; Edmunds 


1996; 


Dwek 


1998: 


Tielens 


1998; 


Edmunds 


2001 


Morgan & Edmunds 


200 


^. In particular 


Eales & Edmunds 


1996). 


Morgan & Edmundsl 


2003 


), and 


Maiolino et al. 


(2004) recognised 



the importance of supernova for the production of the massive amounts of dust observed in 
early galaxies. 



In this paper we focus on the evolution of dust in galaxies that are so youthful that 
AGB stars could not have contributed to the reservoir of dust in these objects. This places 
the burden of dust production solely on supernovae, which in the remnant phase of their 
evolution are also efficient agents of grain destruction. The equations for the evolution of 
dust are therefore greatly simplified, since: (a) they can be formulated using the instanta- 
neous recycling approximation, which assumes that stars return their ejecta back into the 
interstellar medium promptly after their formation; and (b) massive core collapse super- 
novae are the only sources of thermally condensed dust. For the purpose of this analysis we 
have ignored the possible accretion of refractory elements onto preexisting dust in the dense 
phases of the ISM. Consequently, the models developed here allow for the detailed analysis 
of the dependence of the evolution of dust mass on the the dust yield in supernovae and 
the grain destruction efficiency by their remnants in the ISM. Furthermore, since SN play 
the dual role of producing and destroying dust, we also examine the dependence of dust 
evolution on the star formation rate. 

The paper is organized as follows. We first present the analytical models for the evolution 
of interstellar dust in high-redshift galaxies, assuming instantaneous recycling and neglecting 
the contribution of low-mass star to their chemical evolution. We also present analytical 
expressions for the efficiency of grain destruction by SNR in the ISM (§2). The general 
properties of the equations, the evolution of gas mass, galactic metallicity, dust mass, the 
dust-to-gas mass ratio, and their dependence on supernova yields and grain destruction 
efficiencies are discussed in detail in §3. In §4 we apply the general results to the quasar 
J1148+5251 at redshift z = 6.4, when the universe was a mere ~ 900 Myr old. We first 
summarize observations of this quasar, derive its far-IR luminosity and dust mass, discuss 
possible scenarios for the star formation history of this galaxy, and analyze its spectral energy 
distribution (SED), separating it into a starburst, an AGN, and dust emission components. 
The results of our paper are summarized in §5. 
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2. THE EVOLUTION OF GAS AND DUST AT HIGH REDSHIFTS 

2.1. General Considerations 

We define the stellar initial mass function (IMF), 0(m), so that (j){m)dm, is the number 
of stars with masses between m and m + dm, and normalize it to unity in the [mi, m„] mass 
interval, 

rruu 

/ 0(m) dm = 1 (1) 

J mi 

where mi, and m„ are, respectively, the lower and upper mass limits of the IMF. We define 
the IMF-averaged stellar mass, (m), average mass of metals, Y^, and dust. Yd, returned by 
massive stars to the ISM as: 



(m) = / m (j){m) dm (2) 

J mi 

y; = / Y^{m) 0(m) dm/fsN 

J m^D 
^ rmsn 

Yd = Yd{m) (j){m) dm/ f^^ 



where m^ and respectively, the lower and upper mass-limits of stars that becom e 



core-collapse Type II supernovae (SN), taken to be equal to 8 and 40 Mq ( iHeger et al.ll2003l ). 
Yz{m) is the total mass of metals in the gas and dust that is returned by a star of mass m 
to the ISM, Yd{m) is the mass of dust produced by a star of mass m, and /sn is the fraction 
of stars that become SN II, given by: 

rmsn 

/sN = / (t^ijn) dm < 1 (3) 

We define the star formation rate (SFR), ijj{t)j to be the mass of stars formed per unit 
time, and the stellar birthrate, B{t), as the number of stars formed per unit time. For a 
constant IMF the relation between the two is given by: 

m = B{t) (m) (4) 

The Type II supernova rate is given by: 

i?SN(t) = B{t) 0(m) rfm = ^ /sN = ^ (5) 
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where 



JSN 

is the mass of all stars born per SN event. 

In all our calculations we will use three different functional forms for the IMF, charac- 
terized by an m~° power law in the [m;,m„] mass interval. Table [T] lists the parameters of 
the different IMFs considered in this paper. The first functional form is the Salpeter IMF, 
and the choice of the two additional functional forms is motivated by the possibility that 
at high redshifts the stellar IMF may have been more heavily weighted towards high mass 
stars. In addition, the table lists the IMF-averaged values of some relevant quantities. In 
particular, the value of is 147 for a Salpeter IMF, dropping down to 50 Mq for a top 
heavy IMF which has a higher fraction of SN events. 



2.2. The Evolution of the Gas 

Let Mg be the mass of gas in a galaxy. For sufficiently young galaxies we can assume 
that stellar ejecta are instantaneously recycled back into the ISM. In this approximation, 
the rate of change in the ISM mass due to astration and infall and outflow of gas is given 
by: 

MinA [ t \ [Mout 



= -(1 + i::!^ X exp — - X exp 

V Tinf J \ Tinf J \ Tout J \ Tout 

where ip{t) is the star formation rate, R is the fraction of the stellar mass that is returned 
to the ISM, either by SN explosions or quiescent stellar winds. The last two terms are the 
net increase in the mass of the ISM to infall and outflow, respectively. We assumed them to 
be exponential with e-folding times of r^f and Tout, normalized so that the total mass of the 
infalling (outflowing) gas is given by Mj„j (Mout)- 

We assume that the star formation rate is proportional to the gas mass to some power k: 



where Mq is the total mass of the system (stars + gas) at some arbitrary time to- 



-6- 



We will consider two scenarios for the chemical evolution of the system: (1) a closed 
box model in which a galaxy does not exchange any mass with its surrounding medium, that 
is, Minf = Mout = 0; and (2) an infall model, in which a galaxy's mass is initially zero, and 
increases with time due to the infall of metal-free gas. 



2.2.1. Closed Box Model 



In this model the equation for dMg/ dt can then be written as: 



dt 



^0 
Mn 



(9) 



where Mq = Mg{tQ = 0), is the total mass of the system (which initially consists of only 
gas), ipo = ip{t = 0), and Hg{t) = Mg{t)/Mo is the gas mass fraction at time t. Analytical 
solutions are available for arbitrary values of k: 



ljLg(t) = exp 

l-(l-R)(l-k) 



Mn 



k^\ 



(10) 



with the initial condition that /Xg = 1 at time t = 0. 



2.2.2. Infall Model 
In this model the equation for dMg/dt is: 



dMg 
"df 



-(1 - + 



M 



inf 



Hnf 



X exp 



Hnf 



(11) 



Since the initial mass of the system is zero, we define Mq to be the total mass of the galaxy 
at some later time to > 0, which, in the absence of outflows, is given by: 



Mn 



M 



inf 



"Tinf 
Minf 



exp 



t 



'inf 



dt 



(12) 





to Y 


1 — exp ^ 






nnfj. 



Analytical solutions are available for A; = 1 and constant values of Mj„/ and Ti„/, or for 
arbitrary values of k provided that {M^nf / Tinf) ocil){t). 
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Defining fJ,g{t) = Mg{t)/Mo, the solution for the evolution of fJ>g{t) for k 
written as: 



Hg{t) = /i <! exp 
where Tq = Mq/iPq, and 



t 



Hnf , 



— exp 



-{1-R) 



t 

To 



H = 



inf 



;i -i?)VoTin/ - Mo 



1 can be 
(13) 

(14) 



1 — exp 



to 



'inf 



[l-R) 



To 



For a given galactic mass, Mq, there exist only a hmited range of values for ■^q that will 
produce a given gas mass fraction, iJ>g{to), at time to, over a wide range of values for Ti„/. 
The range of values for ipo is given by: 



Hfigito)] Mo X Mo 

S ^0 ^0 S 



l-R 



l-R 



(15) 



where x is the solution for [1 — exp(a;)]/a; = /ig(to)- For values of Tj^/ ranging from r^^j -C to 
to Tinf ^ to, i^o ranges from ~ 130 to 270 Mq yr~^. 



2.3. The Evolution of the Dust 

In very young galaxies with ages < 800 Myr, low-mass stars do not have time to evolve 
off the main sequence and enrich the ISM with their dusty ejecta. We will assume that any 
infalling gas is dust-free, and that the dust only condenses in SN II and does not grow by 
accretion of metals onto grains in the ISM. The equation for the evolution of the mass of 
dust, Md in the galaxy is then given by: 

TI = -^d W\t) + ydRsN[t) (16) 

at Td 

where = M^/Mg is the dust-to-gas mass ratio, is the average yield of dust in Type II 
supernovae, Rsn is the supernova rate in the galaxy, and is the lifetime of the dust against 
destruction by SN blast waves. The lifetime for grain destruction is given by (Dwek & Scalo 
1981; McKee 1987): 

^ M,(t) ^ Mg{t) 

{md) Rsn {misu) Rsn 
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where (m^) is the total mass of elements that are locked up in dust and returned by a single 
supernova remnant (S NR) back to the gas phase eit her by thermal sputtering or evaporative 
grain-grain collisions (jjones et al.lll996l : I Joned 120041 ). and (misM) = /-^d is the effective 
ISM mass that is completely cleared of dust by a single SNR, more formally defined in §2.6 
below. 



Using the expression for the SN rate given by eq. (5), the equation for the evolution of 
dust can now be written as: 



dt 



1 + 



("^ism) 



;i8) 



Solutions for Mdit) depend on the assumed chemical evolution model. 



2.3.1. Closed Box Model 



In this model we can change variables from t to and, using the expression for dfig/dt 
[eq. IQ], rewrite equation (|T8|) in the following form: 



dMd 
dfi 



M, 







, Md-Yd- 

ligj (1 - R) 



where 



u = 



{misu) + 



(19) 



(20) 



:i - R) 

The solution to eq. (IT^ gives the evolution of dust mass (or dust mass fraction) as a function 
of the fractional gas mass: 

Mo 



Md{ilg) 



Yd 



.(?^ism) + Rm^ 
Yd 

{misu) + Rm^ 



(1 - /^r^) 



(21) 



where fid = Md/M^, and where we assumed that Md = at time t = 0, when fig = 1. 
Equation fl2Tl) can be rewritten to give the dust yield required to obtain a given dust-to-gas 
mass ratio, Zd, at a given gas mass fraction fig-. 

{misu) + Rm^ 



Yd 



1 - 



(22) 



A brief glance at eq. fl2Tl) might suggest that the dust mass scales linearly with Mq, 
the total mass of the galaxy. However, the dependence of Md on Mq is more complex, and 
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depends on the amount of grain destruction in the galaxy. When grain destruction is efficient 
so that (misM)^ R rn^, the values of z/ <^ 1 and fi^^^ ^ 1, and eq. fl2Tl) approaches the 
asymptotic value of: 

^ Y.. ^ (23) 

which is independent on the initial mass of the galaxy. The dust yield required to obtain a 
given dust-to-gas mass ratio then becomes: 

Yd = (misM) (24) 

As we will see in §4.3, these asymptotic limits are important when the total mass of the 
galaxy is very uncertain, but dust and gas masses are fairly well determined. 



2.3.2. Inf all Model 

In the infall model, the solution for Mdit), or equivalently /x^, cannot be simply expressed 
as a function of fJ^git). It has an explicit time dependence and is given by: 



l^d{t) = I — 1 



A exp 



'inf 



To 



B exp R)—] +{B - A) exp — 



To 



(25) 



where /i is given by eq. ([H]), tq = Mq/i/jq, and 



^ = (migM)/m^ + l 

B ^ [i-{l-R)]-' 



(26) 



2.4. The Evolution of the Metallicity 

The evolution of the mass of metals, Mz{t), (elements heavier than helium) can be for- 
mally obtained by substituting Yz instead of Yd and by letting (misM) — > in equations (l20l) 
and ([21]) for v and Mdit). 



2.4.1. Closed Box Model 
The solution for Mz{t), and fiz{^^g) = Mz{t)/Mo in this model is given by: 
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Mn 



(1 - ^r^) 



(27) 



where 



;i-i?) 



(2^ 



The metalhcity of the gas, Zg, defined as the mass ratio of metals and the interstellar 
gas is then given by: 



R 



(29) 



2.4.2. Inf all Model 

Substituting into Y^ and setting (misn) to in the solutions of Ma{t) in the infall 
model gives the following equation for fj,z(t) = M^(t)/Mo: 



Y,. 



m^, 



A' exp ( —] - B' exp ( -(1 - R)-] + {B' - A) exp f 



where A' = Tinf/innf - tq), and B' = R ^ 



(30) 



The metallicity of the gas, Zg{t), is given by the ratio: ^zit)/ fig(t), where /ig(t) is given 
by eq. (DSD. 



2.5. The Dust to Metal Mass ratio 

An important quantity is the fraction of the mass of metals in the ISM that is locked 
up in dust, fd, defined as: 



fd = 



Md 



(31) 



- 11 - 



For the closed box model, fd is given by: 



R 



(?^ism) + Rm^ 



,u-l 



(32) 



For the infall model, fd is given by the ratio fid(t)/ fiz(t), given by eqs. ( l25ll and ( l30ll . respec- 
tively. 

An upper limit to this quantity can be obtained by assuming that the dust mass is not 
altered in the ISM, either by grain destruction in SNRs or by accretion in molecular clouds. 
The value of fd is then identical in the infall or closed box model, and given by: 

fd = ^ no grain destruction or accretion in the ISM (33) 

which is the mass fraction of metals in the SN ejecta that condensed into dust grains. Table [2] 
gives the values of F^, Yd for different dust compositions, the total dust yield Y^, and the mass 
fraction of metals that form dust in SN ejecta as a function of stellar m ass for two different 



i nitial metallicities of the progenitor stars. Stellar yields were taken from lWoosley fc Weaver 



(119951 ). The results are also shown in Figured] for an initial stellar metallicity of Zg = Zq. 
The calculations assume that all refractory elements form dust grains with a condensation 
efficiency of unity. To estimate the amount of oxygen that gets incorporated into the dust 
we assumed that the silicate dust grains are of the form [Mg, Fe]2Si04, and the Ca-Ti-Al 
oxides are of the form: Si02, Ti02, and CaO. 

The results show that the mass fraction of metals that can form dust actually decreases 
with stellar mass, as more massive stars form larger amounts of oxygen which is less efficiently 
incorporated into dust. The IMF-averaged mass fraction of metals that can be locked up in 
dust is less than ~ 0.40 (see Table 1). This is a strong upper limit, since these results ignore 
the effect of grain destruction. An observed dust mass fraction of this magnitude in a galaxy 
in which grain destruction is important must therefore imply that the dust mass must be 
enhanced by additional processes, such as the accretion of metals onto grains in the ISM. 



2.6. The Lifetime of Interstellar Dust Grains 

2.6.1. Homogeneous Interstellar Medium 

An important input parameter governing the evolution of the dust is its lifetime against 
destruction by SNR. The mass of ISM gas which has been cleared of dust by a single SNR 
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Fig. 1. — The yield of metals and dust in massive stars for Z — Zq, using the yields of IWooslev fc Weaver 
(|l995l ). Dust yields assume that all refractory elements precipitate out of the ejecta with a condensation 
efficiency of unity. 
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is a measure of the efficiency of this process, and is given by: 

dM 



{misu) = I Cdivs 

VQ 



dv.^ 



dvs (34) 



where Cd{i^s) is the fraction of the mass of dust that is destroyed in an encounter with a 
shock wave with a velocity Vs, {dM/dvs)dVs is the ISM mass that is swept up by shocks in 
the [vs, Vs + dvs] velocity range, and vq and Vf are the initial and final velocities of the SNR. 

Chioffi, McKee, & Bertschinger (1988; hereafter CMB88) presented analytical solutions 
for the evolution of a SNR, expanding into a uniform ISM. Initially, the remnant expands 
adiabatically and its evolution is described by the Sedov- Taylor solution. When the cooling 
time of the shocked gas becomes comparable to its dynamical timescale, the remnant evolves 
as a pressure-driven snowplow (PDS). For sake of our analysis we cast the standard solutions 
in the form of the mass of the swept-up ISM as a function of shock velocity. The evolution 
of the remnant mass, Mgnr, in a medium with solar metallicity can then be written as: 

M,„,(M0) = 400 ^5^6 ^0 0-28 t;;" (35) 

-2.0 for v^ > 1.0 
-1.28 for < 1.0 



a 



dvs 

where £"51 is the energy of the explosion in units of 10^^ erg, no the density of the ISM in 
cm~^, and f^, = fs/wpDS; where fpos is the velocity of the remnant when it transitions from 
the Sedov- Taylor to the PDS phase of its evolution. For a gas of solar abundances and a 
He/H number ratio of 0.10, fpos is given by (CMB88): 

WPDs = 413 nj'^ eI{^^ km s"^ (36) 



Figure [2] depicts the velocity dependence of Cdivs) and the value of (misivi) as a function 
of ISM density. For values of ^ 0.1 — 1.0 cm~^, corresponding to the average density of 
the Galactic ISM, we get that (misM) ~ 1100 — 1300 Mq for an equal mix of silicate and 
graphite grains. For a Galactic ISM mass of Mg = 5 x 10^ Mq, and SN rate of 0.01 yr~^, we 
get an average dust lifetim e ~ 4 x 10^ yr, iri good agreement w ith the value of ~ 4 x 10* yr 
(f« 6 X 10* yr) derived by iJones et al.l (119961 ) and iJoned (120041 ) for the average lifetime of 
silicate (carbon) dust in the Milky Way. 
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Fig. 2. — Left panel - the fraction of the mass of dust that returned to the gas phase by after being 
swept up by a shock of velocity Vg as a function of shock velocity [after Jones et al. ( 1996h ]. Right panel - 



the mass of the ISM gas that is cleared of dust, (misM), is plotted against the density of the homogeneous 
ISM into which the SNR is expanding. The figure depicts the value of (migM) for silicate (dotted line) and 
carbon (solid line) grains. The relation between grain lifetime and (misM) is given in eq. (fT7|) . 

2.6.2. Inhomogeneous Interstellar Medium 

The ISM in galaxies is inhomogeneous, and when the SFR is sufficiently high, as the case 
is for J1148+5251, it is dominated by a hot and low-density gas created by the expanding 
SN blast waves. The multi-phase ISM is then characterized by hot {h), warm {w), and cold 
(c) phases with volume ffiling factors and densities of fi and pi, with ^ /j = 1 (i={h,w,c}). 
SN blast waves propagate predominantly through the low density intercloud medium. The 
shocked phases are in rough pressure equilibrium so that the velocity of the shock propagating 
through, say, the warm phase is related to its velocity in the hot phase by: phV^ ~ Pw^lj. 
The value of (^«ism) for a 3-phase ISM can be written as: 



{misM) = ^ fi fdivjxi) 



dM 



dv^ 



dvs (37) 

-{h,w,c} 

where x is the density contrast between the phase "i" and the dominant hot phase into which 
the remnant is expanding. Since the density contrast between the warm or cold phases to 
that of the hot ISM can be quite large, the shocks propagating through these clouds are quite 
ineffective in destroying the dust (Jones 2004). For a density of ri/i ^ 3 x 10"'^ cm~'^ and a 
density contrast x ^ 10'^(> 10^) between the warm (cold) and hot phases, grain destruction 
in these phases is neglible. 



However, the warm and cold phases are ultimately cycled through the hot phase of the 
ISM by cloud evaporation, cloud crushing by shocks, or cloud disruption by star formation. 
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Injected into a hot 10^ K' 
sputtering on a time scale of (IDwek et al.lll996l : lJonesll2004l ): 



ISM, a dust gra i n of radius a will be destroyed by thermal 



At 



sput 



10' 



a(/im) 



yr 



(3^ 



A grain of radius a = 0.1 fim will therefore survive for a period of Atgput ~ 3 x 10^ yr. In 
high redshift galaxies with total gas masses > 10^° Mq, and star formation rates in excess of 
10^ Mq jt~^, the time scale for the disruption of the cold molecular clouds by star formation, 
~ Mg/ip > 10'' yr, comparable to /Sigput- The effective lifetime of the dust in the 3-phase 
ISM of these objects is therefore determined by Atsput- Using eq. (ITTjl this lifetime can be 
expressed in terms of (misM), giving (misivi) ~ 50 Mq, for Mg = 10^^ Mq, ip = 10^ Mq yr~-^, 
and m^= 150 Mq. 

To keep the model results most general, we will consider the grain destruction efficiency 
as an unknown, and adopt (misn) as a free parameter of the model ranging from (misivi) = 
(no grain destruction) to a value of (misivi) = 1000 Mq. 



3. GENERAL RESULTS 

Figures 3-6 depict the evolution of various quantities as a function of fractional gas 
mass, fig, for closed-box and infall models as a function of time. Results are presented 
for different values of (w^ism)? ranging from to 10^ Mq, corresponding to the range of 
uncertainty in the lifetime of the interstellar dust grains. The value of R, the IMF-averaged 
fraction of the stellar mass that is returned to the ISM over the stellar lifetime, is taken to 
be 0.50. Two different functional forms were used for the IMF: a Salpeter and a top-heavy 
IMF. The IMF parameters and the values of relevant IMF-averaged quantities are given in 
Table [TJ The values of the IMF-averaged dust and gas yields, and respectively, are 
given in Table [H 



3.1. The Evolution of the Gas 

Figure [3] depicts the evolution of the gas mass fraction as a function of time for the 
closed box and infall models. In the closed box model (left panel), the initial gas mass 
fraction is equal to 1, decreasing as the ISM gas is converted into stars. The calculations 
were performed for a SFR law ^{t) = 'tpo{Mg{t) /Mof with Mq = 5 x 10^° Mq, and k = 1.5. 
The curves are labeled by the value of ipito) at time to = 400 Myr, which is our adopted age 
of J1148+5251 (see §4.1). 
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In the infall model (right panel), described by eq. f|T3l) . the initial gas mass fraction is 
zero. It first increases in time as the galaxy accretes mass from its surrounding, but decreases 
later on, when star formation consumes gas at a higher rate than its rate of replenishment by 
infall. The curves are labeled by iJjq, which is related to the current SFR by ipQ = fJ^gito) ■0(to)- 
For each value of ipo, eq. (fT3|l was solved for the value of Tin/ that produced the adopted 
values of Mq and fig at the epoch of to = 400 Myr. 




time(yr) time (Gyr) 

Fig. 3. — The evolution of the gas mass fraction, fJ.g{t), as a a function of time for the closed box (left panel) 
and the infall model (right panel). The various curves are labeled by the current SFR, ip{tQ = 400 Myr) 
for the closed box model, and by ipQ for the infall model. Both quantities are given in units of Mq yr"^. 
The dashed lines depict the values of the adopted mass fraction, fig ~ 0.60, and the age of the quasar, 
t = 400 Myr, at z = 6.4. The figure is discussed in further detail in §3.1 of the text. 

Figure [3] is important for reconstructing the star formation history of a galaxy from 
current observations of the fractional gas mass and the star formation rate. The figure is 
used in §4.5 to construct possible star formation scenarios for J1148+5251. 



3.2. The Evolution of the Dust and Metals 

Figure m depicts the evolution of the dust mass, Md{fig) and mass of metals, M^, nor- 
malized to the initial mass, Mq, for various values of (^ism), the mass of ISM gas that is 
cleared of dust by a single SNR as a function of fig. Initially, fig= 1, but decreases as the gas 
is converted into stars. Calculations are presented for two different stellar IMF: a Salpeter 
IMF (left panel), and a top-heavy IMF (right panel). Initially Md = and rises as the 
ISM is enriched by SN- produced dust. However, eventually the gas and dust in the ISM are 
incorporated into stars, and the mass of interstellar dust decreases. 
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Fig. 4. — The evolution of the mass of dust and metals, Md and respectively, both normalized to the 
initial gas mass, Mq, as a function of the fraction of the ISM gas, fig, in a closed-box model for the chemical 
evolution. Calculations are presented for a Salpeter IMF (left panel) and a top-heavy IMF (right panel). 
Curves are labeled by (misM), the mass of ISM gas that is cleared of dust by a single SNR in units of Mq. 
A value of (misivi) = corresponds to no grain destruction. The horizontal dashed line corresponds to the 
Md/Mo value of 0.004 adopted for J1148-I-5251, and the vertical line to the value of /ig at the epoch of 
400 Myr (see Table [5]). 

Figure presents the same quantities for the infall modeL Both figures show the maxi- 
mum values of fid and fi^ attainable with each IMF. Larger values of fi^ and fiz are obtained 
with a top-heavy IMF. The figures also show that without any grain destruction, the mass 
of dust is simply proportional to the mass of metals, but decreases more rapidly than the 
mass of metals when grain destruction is taken into account. 



3.3. The Evolution of the Dust-to-Metals and Dust-to-Gas Mass Ratios 

The effect of grain destruction is to decrease the fraction of condensable elements in 
the solid phase of the ISM. This point is illustrated in Figure |6] which depicts the evolution 
of the mass fraction of metals locked up in dust, fd = Md/M^, versus fig for the closed 
box model. When (misivi) = 0, fd is constant and equal to Yd/Y^ ~ 0.35, the fraction of 
the metals in the SN ejecta that condensed and formed dust. As the figure illustrates, this 
fraction deceases with fig as the grain destruction efficiency increases. Similar quantitative 
results can be obtained for the infall model. 

Figure [7] shows that when grain destruction is ignored, i.e. (misM) = 0, the dust-to-gas 
mass ratio, Zd, continues to rise since both the gas and dust are incorporated into stars. 



- 18 - 




Fig. 5. — The evolution of the mass of dust and metals, Md and respectively, both normalized to the 
gas mass, Mq, at to = 400 Myr as a function of time for the infall model. Calculations are presented for 
a Salpeter IMF (left panel) and a top-heavy IMF (right panel). Curves are labeled by (toism) hi units of 
Mq. The horizontal dashed line corresponds to the Md/Mo value of 0.004 adopted for J1148+5251, and the 
vertical line corresponds to the adopted galaxy's age of 400 Myr (see Table [5]). 

but the ISM is continuously enriched by dust formed in SN ejecta. When grain destruction 
is taken into account, reaches a steady-state at values of fig which become increasingly 
smaller as the grain destruction efficiency, which is related to the value of (misM), increases. 



3.4. The SN Dust Yields Needed to Produce an Observed Dust-to-Gas Mass 

Ratio 

Figure [8] shows how much dust an average SN must produce in order to give rise to a 
given dust-to-gas mass ratio, for various grain destruction efficiencies. The value of 1^ was 
calculated when Hg reaches a value of 0.60, the adopted gas mass fraction of J1148-I-5251 at 
400 Myr. The figure shows that, for example, to produce a value of Zd = 0.0067 at fig = 0.60, 
a SN must produce about 0.4 (1.2) Mq of dust for a top-heavy (Salpeter) IMF, provided the 
dust is not destroyed in the ISM. Even with modest amount of grain destruction, (misn) = 
100 Mq, the required SN dust yield is dramatically increased to about 1 — 2 Mq, depending 
on the IMF. The horizontal line in the figure corresponds to a valu e of = 0.02 Mq, th e 



largest amount of dust directly observed in the ejecta of a supernova (ISugerman et al.ll2006l ). 
The figure shows that even without grain destruction, the largest observed yield can only 
give rise to a gas-to-dust mass ratio of ~ 3 x 10""^. 
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Fig. 6. — The evolution of the mass fraction of metals that is locked up in dust, fd given by eq. (20), is 
plotted versus fig, the fraction of the ISM gas for the closed box model. Curves are labeled by (migM) in 
units of A/q. In the absence of grain destruction the value of fd is equal to ~ 0.35, the fraction of metals in 
the SN ejecta that condenses into dust. 

4. APPLICATION TO J1148+5251 

The results of our chemical evolution model can be readily applied to any galaxy suffi- 
ciently young so that AGB stars are only minor contributors to the dust abundance in the 
ISM. Here we concentrate on the quasar J1148+5251. 



4.1. Observational Properties 

Table [3] summarizes the observed properties of J1148+5251. At redshift z = 6.4 the age 
of the universe is 890 Myr for a ACDM universe with Q^n = 0.27, = 0.73 and a Hubble 
constant Hq = 70 km s^^ Mpc^^. 

Figure [9] depicts the observed far-IR and submillimeter fluxes at the observed wave- 
lengths. The different curves are spectral fits to these fluxes for different dust compositions. 



The o ptical properties for the silicate and graphite grains were taken from iDraine &: Lee 



(Il984j ) and for the carbon grains from Rouleau & Martin (1991). Dust masses vary from 
~ 10^ to 5 X 10* Mq, depending on dust composition. Table H] summarizes the derived 
properties and IR luminosities for the different dust compositions. 

For comparison with model calculations we assumed that the onset of star formation 
in this galaxy occurred at z = 10, when the universe was 490 Myr old, giving an age of 
400 Myr for this object at z = 6.4. Based on the CO observations (see Table [3]) we adopted 
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Fig. 7. — The evolution of the dust-to-gas mass ratio, Z^^ and the dust-to-metals mass ratio, as a function 
of [ig. Curves are labeled by (misM) in units of Mq. The dashed horizontal and vertical lines represent, 
respectively, the adopted gas-to-dust mass ratio and gas mass fraction of J1148-f 5251 at t = 400 Myr. The 
figure shows that when grain destruction is important, (?nisM)?J 100 Mq, a top-heavy IMF is required to 
produce the observed amount of dust at p^g = 0.6. 

a molecular gas mass of 1.5 x 10^° M©, and and an equal mass of atomic gas, for a total 
gas mass of 3 x 10^° M©. From the far-IR emission we adopted an average dust mass of 
Mrf = 2 X 10^ Mq. From estimates of the dynamical mass, we adopted a total galactic mass 
of Mo = 5 X 10^'' Mq at time t = to = 400 Myr. This mass corresponds to the initial mass 
of the galaxy in the closed box model, and to the total mass accreted until to in the infall 
model. With these values, we get a gas mass fraction of fig = 0.60, and a dust-to-gas mass 
ratio of Za = 0.0067 at t = 400 Myr. Finally, we adopted a far-IR luminosity of 2 x 10^^ Lq. 
Table [5] summarizes the various derived and adopted properties of J1148-I-5251. 



4.2. The Dust Mass and the Required Dust Yield in Core-Collapse SN 

Figures m and Ogive an upper limit on the mass of dust that can be produced in 400 Myr 
by supernovae when grain destruction is ignored, that is, when (misM) = 0. The maximum 
dust mass is ~ 10^ Mq (/id 2 x lO'^) for a Salpeter IMF, and ~ 5 x 10^ M© (/i^ ^ lO'^) for 
a top-heavy IMF. The observed dust mass of 2 x 10^ Mq requires therefore a top-heavy IMF, 
which can produce this mass of dust even with a grain destruction efficiency corresponding 
to a value of (misM) = 100 Mq. This value of (rriisM) corresponds to a grain lifetime of: 

r. = , = 1^ = 27 My. (39) 

(misM) iisN 100 X 7.4 ' ^ ' 
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Fig. 8. — The IMF-averaged yield of dust by type II supernova, Y^, that is required to account for a given 
dust-to-gas mass ratio Z^,, is presented for different values of (misM) given in units of Mq. The horizontal 
dashed line near the bottom of the figure corresponds a value of 1^ = 0.02 Mq, the highest inferred yield 
of dust in a supernova to date (Sugerman et al. 2006). The vertical dotted line represents the value of 
at = 0.60. Curves are labeled by (misM) given in units of Mq. Solid and dashed lines correspond to 
calculations done for a top-heavy and a Salpeter IMF, respectively. The top two dashed (solid) horizontal 
lines represent IMF-averaged theoretical dust yields for a Salpeter (top-heavy) IMF. The results are identical 
for both, the closed box and infall models. 
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Fig. 9. — Spectral fits of several possible dust compositions to the observed far-IR and submillimeter 
observations of J1148+5251. References to the observations are listed in Table [3l Temperatures, masses, 
and luminosities of the different dust types are listed in Table HI 

where the SN rate was calculated for a SFR of 10'^ Mq yr~^and a Salpeter IMF for which 
(msN) = 147 Mq. 

This result is also illustrated in figure [8] that shows that the yield of dust per SN required 
to produce a dust-to-gas mass ratio of 0.0067 with (misM) = 100 Mq is about 1 Mq per SN. 
This is about equal to the total mass of condensable elements produced in a SN with a 25 Mq 
progenitor star (Woosley & Weaver 1995; Nomoto et al. 2006). Such average dust yield is 
theoretically only attainable with a top-heavy IMF, provided that the dust condensation 
efficiency in the SN ejecta is about 100%. Observationally, this yield is significantly higher 
than the ~ 0.02 Mq of dust in SNII 2003gd in the galaxy NGC 628, i nferred from the 



analy sis of the IR emission and the internal extinction in the SN ejecta ( jSugerman et al. 



20061 ). If this low yield is typical, then SNe cannot be the dominant source of dust in these 



young galaxies. Alternative mechanisms that can produce the observed dust m ass, such as 



accretion in molecular clouds, or formation around the AGN (lElvis et al.ll2002l ). need then 
to be included in the dust evolution model. 
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4.3. The Dependence of the Required SN Dust Production Yield on the Total 

Mass of the Galaxy 

The dust yield required to produce a given dust-to-gas mass ratio at time to depends 
on the total mass of the galaxy, which we took to be the dynamical mass of J1148+5251. 
Since this quantity is uncertain, we explore the dependence of the dust mass produced in 
the model on the adopted total mass of the system. For the closed box model, the relation 
between these two quantities, and Mq, is given by eq. (!2T|) . 

The total mass of the galaxy must exceed its gas mass, which we took to be 3 x 10^'' Mq, 
so we varied Mq from 4 to 8 x 10^° Mq, which spans the dispersion in the observed dynamical 
mass of the galaxy. Figure [10] depicts the results of our calculations for the two different 
IMFs and the three different values of (misn) used in figure [H In the absence of grain 
destruction Yd decreases by a factor of ~ 4 — 5 over the entire mass range of Mq. However, 
when grain destruction is taken into account. Yd becomes increasingy independent of the 
total galactic mass, approaching the asymptotic behavior for (misM) ^ R nT-* presented in 
eq. (12^ . For example, for (misM) ~ 300 Mq, and a dust-to-gas mass ratio of Zd ~ 0.003 we 
get that Yd> 1 Mq, significantly larger than any dust mass observed in any SN ejecta. 



4.4. The Star Formation Rate in J1148+5251 

From the observed luminosity and gas mass we can derive the current SF R using the 



relatio ns between the SFR and the far-IR luminosity and the total gas mass (IKennicutt 



1998al jb[). For a Salpeter IMF the SFR is related to the far-IR luminosity of the galaxy by: 
^(Mq yr-^) = 1.7 X 10"^° L{Lq) Salpeter IMF (40) 
which for a value of Ljr = 2 x 10^^ Lr^ gives a SFR o f ~ 3 400 Mq yr^^ A high SFR 



of ~ 3000 Mq yr ^ was also derived by iMaiolino et al.l (120051 ) from the luminosity in the 
[C II] 158 fim line detected in this quasar. 

The SFR can also be derived from the gas mass using the empirical relation between 
the SFR per unit area, Ssfr, and the gas mass surface density, in star forming galaxies 
(Kennicutt et al. 2006): 

/ y, N 1.56±0.04 

EsFR (Mq yr-1 kpc-^) = {^.b^l'l) x 10"^ ' , Salpeter IMF (41) 

\M0 pc / 

The CO maps of J1148+5251 suggest that it consists of two blobs of comparable size with 
a diameter of about 1 kpc. For an adopted gas mass of 3 x 10^° Mq we get a total SFR 
ip ^ 300 Mq jt'^. 
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Fig. 10. — The dependence of the dust yield in an average Type II supernova required to produce the 
observed dust mass of 2 x 10^ Mq, on AIq, the total mass of the system (see eg. [^Tt . The calculations were 
performed for a gas mass of Alg = 3 x lO^*' Mq and a dust mass of Md = 2 x 10^ Mq. Curves are labeled 
by (misM) given in units of Mq. Solid and dashed lines correspond to calculations done for a top-heavy and 
a Salpeter IMF, respectively. When grain destruction is taken into account, te required dust yield becomes 
independent of the total mass of the galaxy. 
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The two different tracers give very different star formation rates. Both estimates are 
highly uncertain: the first assumes that the dust reradiates all the starburst's luminosity 
and ignores any possible contribution of an AGN to the heating of the dust. The second 
assumes that the CO observations traces the densest gas regions in the object, as inferred 
from HCN-CO correlat ions, and is therefore a good measure of the total gas mass in galaxies 
( iGao fc Solomon! I2OO4I ). However, HCO"'" observations which probe denser regions of molec- 
ular gas, have cast doubts on the r eliability of HON (and hen ce CO) as an unbiased tracer 
of dense molecular gas in ULIRGs (IGracia-Carpio et al.l 120061 ). 



Another major source of uncertainty is the stellar IMF. For example, the SFR derived 
from the IR luminosity using a Salpeter IMF, will decrease from a value of 3400 Mq yr~^, to 
a rate of about 380 Mq yr~^ for a mass- or top-heavy IMF. Because of all these uncertainties, 
the star formation rate of J1148+5251 cannot be uniquely determined. In the following, we 
discuss the dependence of its star formation history on the assumed current SFR. 



4.5. The Star Formation History of J1148+5251 

In the closed box model, the evolution of the dust, and metals could be presented as a 
function of the gas mass fraction iig (Figures HHT]). The translation of the dependence of these 
quantities from /i^ to time requires knowledge of the star formation history of J1148+5251, 
which in turns depends on the assumed SFR and the initial gas mass of the object. For 
example, given an initial gas mass of Mq = 5 x 10^° Mq, a SFR of 3000 M© yr~^ will 
exhaust almost all the available gas in less than ~ 100 Myr. This suggests one or more of 
the following possibilities: (1) the onset of star formation occurred relatively shortly before 
the observations; (2) the galaxy started with an initially larger reservoir of mass; (3) that 
rapid infall, comparable to the SFR, kept the reservoir of gas sufficiently high. 

Some of these possibilities were presented more qualitatively in figure [3l which depicts 
the evolution of the mass fraction of gas as a function of time for the closed box and infall 
models. For the closed box model (left panel), the calculations assume an initial gas mass 
of 5 X 10^° Mq, and the Kennicutt law for the relation between the SFR and the gas mass: 
ijjit) oc Mg(t)^, with k = 1.5 [see eq. ([7])]. The solid line depicts the evolution of figit) for an 
initial SFR ipQ = 150 Mq yr~^, corresponding to a value of ■?/' = 70 Mq yr~^ at t = 400 Myr, 
which reproduces the adopted value of fig = 0.60 at that epoch. The additional curves depict 
the evolution of figit) for initial SFRs of 650, 2.2 x 10^, and 6.5 x 10'^ Mq yr~^, corresponding 
to values of tpito) = 300, 1 x 10^ and 3 x 10^ Mq yr'\ at to = 400 Myr. The corresponding 
gas mass fractions at that epoch are: 0.19, 0.036, and 0.0052, respectively. The figure 
shows that a low observed SFR of only 70 Mq yr~^ is required to fit the observations of 
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J1148+5251, given the initial conditions and assumptions summarized in Table [51 A current 
SFR of 3000 Mq yr~^ requires changes in the initial conditions and model assumptions. If 
the galaxy is indeed 400 Myr old, then the value of Hg is 0.0052, requiring the initial gas 
mass of J1148+5251 to be 5 x 10^°/0.0052 ~ 1 x 10^^ Mq. A more plausible scenario is that 
fig is 0.60, but that the age of the starburst is only about 10"^ yr. 

Similar conclusions are reached for the infall model. The right panel of Figure [3] depicts 
the evolution of the gas mass, constrained to fit the adopted values of Mq and Hgito) at 
time to = 400 Myr. The fit requires the values of the product ifjQ to to be between ~ 
(5 — 11) X 10^° Mq. If star formation has been an ongoing process over a period of 400 Myr, 
then the current SFR must be between ~ 125 and 285 Mq yr~^. To accommodate a much 
larger SFR, say of 3000 Mq yr~\ requires the age of the starburst to be about 3 x 10^ yr. 



4.6. The Spectral Energy Distribution of J1148+5251 

The SED of J1148+5251 offers very few clues regarding the relative starburst or AGN 
contribution to the thermal dust emission from this galaxy. Figure [11] depicts the galaxy's 
SED from UV to submillimeter wavelengths. The UV and optical parts of the spectrum are 
most likely dominated by escaping starlight, and the far-IR by reradiated thermal emission 
from dust. The dashed-dotted grey l ine in the figure represents the intrinsic stellar radia- 



tion field synthesized with PEGASE ( iFioc fc Rocca-Volmerangelll997l ) for a continuous star 
formation rate of age t = 400 Myr, with a top-heavy IMF, and a SFR of 2500 Mq yr"^. 
The total intrinsic stellar luminosity is 1.0 x 10^^ Lq. Part of this stellar energy is absorbed 
by dust and re radiated at IR wa vlengths. We used a simple screen model with a Galactic 



extinction law (IZubko et al.ll2004j ) to calculate the spectrum of the escaping stellar radiation, 
depicted by the dotted line in the figure. The magnitude of the extinction was chosen so that 
the total energy absorbed by the dust, shown as a grey shaded area in the figure, is equal to 
the total reradiated far-IR emission. The total luminosity radiated by the starburst-heated 
dust is 4.6 X 10^^ Lq. The composition of this dust was taken to consists of a mixture 
of silicate and graphite dust with mass fractions of 2/3 and 1/3, respectively, and a 
distribution of dust temperatures ranging from 40 to 150 K. 

The thick solid line represents the sum of all emission components, and represents the 
best fit of select model parameters (the intensity of the starburst, and the slope of the 
power law describing the AGN spectrum) to the observations. The model described above 
is only a plausible one, and definitely not unique. The possible existence of many distinct 
emission components, and the uncertainty in the IMF illustrate the difficulty in determining 
the star formation rate from the galaxy's SED. 
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The origin of the near- to mid-IR (NMIR) emission is more uncertain. The rest-frame 
~ 0.5 — 1 yum fluxes are in excess of the stellar emission that can be produced by a young 
starburst. It also cannot be produced by dust, since it requires the grains to radiate at 
temperatures above their sublimation point of ~ 1500 K. We therefore fit the near- to mid- 
IR emission with two components: an AGN represented by a i>~^'^ power law, and a hot 
dust component represented by a 1:9 mix (by weight) of silicate and graphite grains with 
a distribution of dust temperatures ranging from 150 to 1500 K. The relatively low 
silicate-to-graphite mass ratio was chosen to avoid the production of a mid-IR excess due to 
the 9.7 fim silicate emission feature. The power law is in good agr eement with the a verage 
value found in the sample of ISO Palomar-green QSOs studied by iHaas et al.l (120031 ). The 
total intrinsic luminosity of the AGN is about 7 x 10^'^ Lq, and the luminosity radiated by 
the AGN-heated dust is 1.3 x lO^^ Lq. 



If the black hole (BH) radiates at the Eddington luminosity: 

'Mbh 



-leddX 



3 X 10" 







(42) 



then the BH mass required to pr oduce the AG N lumi nosity is ~ 2 x 10^ M©, comparable 
to the mass estimate derived by IWillott et al.l (120031 ) from the width of the Mg II line. 
Mechanisms for the formation of seed black h oles that enable their rapid growth to masses 
in excess of ~ 10^ M© have been discussed by lLodato fc NatarajanI (120061 ). 



5. SUMMARY AND DISCUSSION 

The early universe is a unique environment for studying the role of massive stars in the 
formation and destruction of dust. In this paper we developed analytical models describing 
the evolution of the gas, dust, and metallicity in high redshift galaxies. The equations 
describing their chemical evolution can be greatly simplified by using the instantaneous 
recycling approximation, and by neglecting the delayed contribution of low mass stars to the 
metal and dust abundance of the ISM. Neglecting any accretion of metals onto dust in the 
interstellar medium, the evolution of the dust is then completely driven by the condensation 
of refractory elements in the ejecta of Type II supernovae, and the destruction by SN blast 
waves in the interstellar medium. The solutions for the evolution of the mass of gas, dust, 
and metals are presented in §2 for closed box and infall models for the chemical evolution 
of the galaxy, and for different functional forms for the stellar initial mass function. The 
results of our paper can be briefly summarized as follows: 

1. The maximum attainable dust-to- metal mass ratio in any system is equal to the IMF- 
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Fig. 11.— SED of J1148+5251. The observed fluxes from J1148+5251 are plotted as a 
function of the rest frame wavelength for a redshift of 6.4. Data and references are given 
in Table [61 The galaxy's SED consists of four distinct emission components: the energy 
released by the starburst; the emission from the AGN; the starlight absorbed by the dust 
and reradiated at far-IR wavelengths; and the emission from the AGN that is absorbed by 
the dust and reradiated at mid-IR wavelengths. The gray area in the flgure represents the 
stellar energy absorbed by the dust, and the hatched area the fraction of the ionizing stellar 
radiation that is absorbed by the gas. Details on the derivation of the emission components 
are given in §3.5 of the text. 
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averaged mass fraction of metals that are refractory and able to condense onto grains 
in SN ejecta, which is about 0.35 (see Table [2]). 

2. This mass fraction is significantly reduced when grain destruction is taken into account 
(Figure [6]). An observed dust-to-metals mass fraction > 0.4 will therefore imply that 
accretion of ices onto interstellar grains in the ISM may be important in determining 
the dust mass in the galaxy. 

3. Grain destruction plays an important role in the evolution of dust. However, its effi- 
ciency depends on the morphology of the ISM and is therefore highly uncertain (§2.6). 
We therefore present all our results for different values of (misM), the effective mass of 
ISM gas that is completely cleared of dust by a single SNR. 

4. In §3 we present the general results of our models, describing the evolution of the gas, 
the dust, and the metals for both, the closed box and infall models. 

5. We applied our general results to J1148+5251, a dusty, hyperluminous quasar at red- 
shift z = 6.4. The observed and adopted quantities of J1148+5251 are summarized in 
Tables [3111 and El 

6. The formation of about 2 x 10^ Mq of dust in this galaxy requires an average SN to 
produce about 1 Mq of dust (Fig. [HI). Theoretically, such large amount of dust can be 
produced if stars are formed with a top-heavy IMF, and with a moderate amount of 
grain destruction ((misM)~ 100 Mq). A Salpeter IMF fails to produce this amount of 
dust even in the absence of any grain destruction. Observationally, the required dust 
yield is in excess of the largest amount of dust (~ 0.02 Mq) observed so far to have 
formed in a SN. This suggests that accretion in the ISM may play an important role 
in the growth of dust mass. 

7. Figure [Til depicts the galaxy's spectral energy distribution from UV to far-IR wave- 
lengths. The SED includes emission from the starburst, the AGN, and hot and cold 
dust components, radiating at mid- and far-IR wavelengths, respectively. 

8. Uncertainties in the fraction of the infrared luminosity that is powered by the starburst 
and in the stellar IMF prevent any accurate determination of the current star formation 
rate in the galaxy, or the unique determination of its star formation history (see §3.4 
and Fig. [3D. 

9. Simple decomposition of the galaxy's SED into its emission components suggest that 
the intrinsic starburst luminosity is about 1 x 10^"^ Lq, 4.6 x 10^^ Lq of which is 
absorbed and reradiated by dust at far-IR wavelengths. 
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10. The ~ 3 //m IR emission from the galaxy can neither be produced by starhght nor 
hot dust. It therefore must be emission from the AGN, and we estimate the AGN 
luminosity to be about 7 x 10^^ L©, 1.3 x 10^^ Lq of which is assumed to be absorbed 
and reradiated by dust at mid-IR wavelengths. 

11. The AGN luminosity requires the formation of a black hole of a mass ^ 2 x 10^ Mq 
at this redshift. 
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Table 1. Values of IMF-averaged Quantities^ 





















IMF 


a 




TTlup 


(rn) 




/sN 


Y, 


Yd 


Salpeter 


2.35 


0.1 


100 


0.35 


147 


0.0024 


1.4(1.7) 


0.5 (0.6) 


mass-heavy 


2.35 


1.0 


100 


3.1 


58 


0.054 


1.4 (1.7) 


0.5 (0.6) 


top-heavy 


1.50 


0.1 


100 


3.2 


50 


0.064 


2.2 (2.7) 


0.7 (0.9) 



^See §2.1 for the definition of aU quantities. Masses and yields are in units of Mq. 
Metallicity and dust yields were calculated for metaUicities of O.OIZq, and Zq (in 
parenthesis). 



Table 2. Maximum Dust Yield in Massive Stars 



Mass Z/Zq Metals Sihcatcs Carbon Ca-Ti-Al Dust Yd/Y^, 



12. 





01 


0.4701 





2848 





0901 





0043 





3792 





8066 




1 





0.6712 





2697 





0815 





0235 





3747 





5583 


13. 





01 


0.6413 





3670 





1090 





0077 





4837 





7542 




1 





0.7869 





3691 





1151 





0090 





4932 





6268 


15. 





01 


0.9487 





4164 





1490 





0110 





5764 





6075 




1 





1.4340 





4759 





1623 





0248 





6630 





4623 


18. 





01 


1.6982 





4414 





1940 





0162 





6516 





3837 




1 





2.1711 





5337 





2493 





0288 





8118 





3739 


20. 





01 


2.4372 





5689 





2080 





0282 





8052 





3304 




1 





3.0787 





7910 





2143 





0287 


1 


0341 





3359 


22. 





01 


2.9129 





8025 





2480 





0333 


1 


0839 





3721 




1 





3.7869 


1 


1356 





2424 





0414 


1 


4193 





3748 


25. 





01 


4.0769 


1 


0311 





2790 





0333 


1 


3435 





3295 




1 





5.0297 


1 


1484 





3234 





0704 


1 


5422 





3066 


30. 





01 


6.0489 


1 


3282 





3150 





0439 


1 


6871 





2789 




1 





7.2639 


2 


0342 





2916 





1159 


2 


4417 





3361 


35. 





01 


8.3389 


1 


9632 





3520 





0494 


2 


3647 





2836 




1 





9.7242 


2 


5560 





3216 





1590 


3 


0366 





3123 


40. 





01 


10.5787 


2 


3926 





3890 





0722 


2 


8538 





2698 




1 





11.8819 


3 


0309 





3696 





1919 


3 


5924 





3023 



Based on 



Woosley fc Weaver 



([1995|) yields of massive stars. 
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Table 3. Observed Properties of the QSO Jl 148+5152 



Observed quantity 



Value 



Reference 



R.A. (J2000) 
Dec. (J2000) 



ll''48'"16!6 
+52°51'50" 
0.2" 



(Fan et al. 2003) 
(Fan et al. 2003) 
(Walter et al. 2004) 



z(Lya) 
z([Mg II]) 
z(CO) 



6.37 ±0.03 
6.403 ±0.005 
6.419 ±0.001 



(White et al. 2003) 
(Iwamuro et al. 2004) 
(Bertoldi et al. 2003b) 



M(CO(3~2)) 
M(CO(7-6), CO(6-5)) 

Mdyn. 

Mbk 



- 1.6 X 10^° Mq 

- 2 X 10^° Mq 

(5.0 ±2.5) X lO^"^ M(r 

3 X lO'' Mq 



(Walter et al. 2004) 
(Bertoldi et al. 2003b) 
(Walter et al. 2004) 
(Willott et al. 2003) 



Table 4. Dust Properties and IR Luminosities^ 





Tdust (K) 


Mdust (Mq) 


Lm (Lq) 


Graphite 


49 


2.7 X 10*^ 


1.9 X 10^3 


Silicate 


47 


4.9 X 10* 


2.0 X 10^3 


Carbon BE 


64 


9.7 X 10^ 


2.4 X 10" 


Carbon AC 


74 


9.3 X 10^ 


2.9 X 10^3 



^Graphite and sili cate optical properties were taken from 
(jPraine fc Ledll984l ). and the optical prope rties of the car- 
bon dust were taken from (jRouleau &: Martin., 1991.) . 
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Table 5. Summary of Derived and Adopted Properties of J 1148+5251 



Quantity Value 

Age 400 Myr 

Initial Mass (closed box model) 5 x 10^° M© 

Gas mass at 400 Myr 3 x 10^° Mg 

IR luminosity 2.0 x lO^^ Lq 

UV-optical luminosity 2.0 x 10^^ Lq 

Dust Mass 2 x 10^ M0 

Star formation ralo W IViiK) M . v 



^At z = 6.4. 
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Table 6. Observed Fluxes From J1148+5251. 



\ 

'^obs 


\ 

■^rest 




References 




(/im) 


i "I y "1 1 r 1 

[mjj) 




0.77 


0.10 


0.0017 ± 0.0002 


fFan et al. 2003) 


0.91 


0.12 


0.0325 ± 0.0033 


fFan et al. 2003) 


1.08 


0.14 


0.0887 ±0.0089 


fFan et al. 2003) 


1.22 


0.16 


0.0796 ± 0.0080 


fFan et al. 2003) 


1.63 


0.22 


0.091 ±0.009 


fWillott et al. 2003) 


2.19 


0.30 


0.103 ±0.010 


fWillott et al. 2003) 


3.6 


0.49 


0.124 ±0.002 


fJiane et al. 2006) 


4.5 


0.61 


0.140 ±0.003 


r Jianer et al. 2006) 


5.8 


0.78 


0.133 ±0.010 


r Jianer et al. 2006) 


8.0 


1.08 


0.241 ±0.016 


fJiane et al. 2006) 


16 


2.16 


0.51 ±0.25 


fCliarmandaris et al. 2004) 


22 


2.96 


0.74 ±0.37 


fCharmandaris et al. 2004) 


24 


3.23 


1.52 ±0.13 


fJianer et al. 2006) 


70 


9.43 


<io 


fJiane et al. 2006) 


350 


47.2 


21.0 ±8.1 


fBeelen et al. 2006) 


450 


60.7 


24.7 ±8.6 


fRobson et al. 2004) 


850 


115 


7.8 ±0.7 


fRobson et al. 2004) 


1200 


162 


5.0 ±1.1 


fBertoldi et al. 2003a) 


3000 


404 


<0.52 


rSertoldi et al. 2003b) 


2.1 X 10^ 


2.9 X 10^ 


0.055 ±0.012 


fCarilli et al. 2004) 



